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Abstract: We report measurements of nonlinear magneto-optical rotation 
(NMOR) for the D2 line of ^^Rb atoms in an antirelaxation-coated vapor 
cell in the presence of a radio-frequency (rf) field. The experimental NMOR 
signals as a function of rf field frequency for various rf field powers are 
compared to a theoretical model based on the density-matrix formalism. 
The comparison between experiment and theory enables understanding of 
the ground-state atomic spin polarization dynamics, illustrated using plots 
of the probability distribution of the atomic angular momentum. 
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1. Introduction 

The dynamics of atomic spin polarization in the presence of radio-frequency (rf) fields has 
been extensively studied (for example, see lH] 121 13] IH). Recently, it has been demonstrated 
that nonlinear magneto-optical rotation (NMOR) is a sensitive probe of atomic spin dynamics 
|l5]|6]|7l, and can in fact be used to selectively create and probe different multipole moments 
of the atomic spin polarization |[8l|9][T0l. There are numerous practical applications for the use 
of NMOR to detect the response of atomic spins to rf fields: for example, in nuclear magnetic 
resonance (NMR) [11 j . including nuclear quadruple resonance (NQR) I.12J , and magnetic reso- 
nance imaging (MRI) ifTSll experiments. One can also use such methods in tests of fundamental 
physics fT4]. In our previous research, we built an alkali-vapor magnetometer for the detection 
of rf fields I.15il . Recent work by researchers using a similar experimental setup is described in 
Refs. |fT6l[T7]| . The experiment ifTSll demonstrated a sensitivity to oscillating magnetic fields of 
100 pG/-\/Hz. The line shapes observed in the magnetometry experiment in the limit of low rf 
power are well understood. However, when higher-strength rf fields are applied, nontrivial line 
shapes are seen, with similarities to those of the of the "Majorana-Brossel effect" in double- 
resonance spectroscopy fTF|. We have conducted experimental and theoretical investigations to 
explain the mechanism that produces these line shapes. 

In this work we present experimental NMOR signals taken on the D2 line of ^^Rb atoms 
contained in an antirelaxation-coated cell in the presence of an rf field. The signals are primarily 
due to interaction with the Fg = 2 ^ Fg = \ transition (the subscripts g and e indicate the 



ground and excited states, respectively). They are compared to a theoretical model based on the 
density-matrix formalism. The theoretical model, which is found to be in good agreement with 
experimental results, enables the understanding of the underlying atomic spin dynamics and 
their relationship to the detected NMOR signals. The character of the observed line shapes is 
found to be different depending on the value of the Rabi frequency for the rf field, Hif, relative 
to the ground-state atomic-polarization relaxation rate 7, and the magnetic resonance frequency 
D.L- This work is motivated by an ongoing project aimed at measuring collisional transfer of 
alignment in collisions between different ground-state alkali atoms ||T9l . 



2. Description of the experiment and theory 

2.1. Description of experhnent 

The experiment employs a spherical paraffin-coated glass vapor cell (diameter = 10 cm) fiUed 
with a natural isotopic mixture of rubidium. The cell coating allows polarization of ground-state 
alkali atoms to survive several thousand wall collisions [20, .21 J. thereby extending the lifetime 
of atomic polarization. The cell is heated by a constant air stream from a heat exchanger. Dur- 
ing all measurements in this work the temperature of the cell is stabilized at « 41°C, resulting 
in a Rb vapor density of 7.2(7) x lO"^ atoms per cm^. The density was extracted from fitting a 
low-light power (« 8 |jW) transmission spectrum for the Rb D2 line to a calculated spectrum 
assuming linear absorption. Under these experimental conditions the longitudinal relaxation 
rate of ground-state alignment of the isotope ^^Rb was measured to be 31.7(7) s^'. The value 
represents the effective relaxation rate of alignment polarization due to diverse types of atomic 
collisions occurring in the vapor cell, e.g. electron-randomization collisions with the wall, uni- 
form relaxation due to the reservoir effect and spin-exchange collisions among Rb atoms Ii22l . 
A four-layer p-metal magnetic shield (see Fig. [TJ surrounds the vapor cell and reduces the ex- 
ternal magnetic field by a factor of ^ 10^ f2V\. Additionally, three orthogonal solenoidal coils 
placed inside the innermost shield are used to compensate residual magnetic fields to below 1 
jjG and to create additional static and oscillating fields. For our measurements we apply a static 
magnetic field Bo in the z direction (typical strength Bq = 0.8 mG) and an oscillating magnetic 
field Brf — Brfxcos Orff along the x-axis (typical amplitude up to Brf — 0.25 mG) using the sine 
wave output of a lock-in amplifier. 

A laser beam, initially polarized along the z-axis, propagates through the vapor cell in the y 
direction. The beam is generated by a distributed feedback (DFB) laser, which is tuned to the 
red wing of the Fg ^2 ^ Ff, ^ 1 transition in the 780 nm D2 line of ^^Rb. The wavelength is 
locked using a dichroic atomic vapor laser lock [23 , 24|. The light intensity is 10 |jW throughout 
the measurements and the beam diameter is w 2 mm. Linear dichroism of the atomic medium 
induces changes in the light polarization. (For high light power, the medium can also acquire 
circular birefringence f25l.) After transmission through the vapor cell the polarization of the 
light beam is analyzed using a balanced polarimeter setup, consisting of a polarizing beam 
spUtter (Rochon crystal) and two photodiodes detecting the intensities of the two beams exiting 
the crystal. The component of the difference signal that oscillates at the frequency of the rf field 
is then extracted by the lock-in amplifier 

Experimental signals as a function of rf frequency, along with predictions of the theory de- 
scribed in Sec. 2.2 are given in Fig. |2] for different amplitudes of the rf field. At the lowest 



amplitudes, the observed fine shapes are Lorentzians, while for higher amphtudes, additional 
features are seen. 



2.2. Description of theory 

The experimental signal is primarily due to interaction with the Fg ^ 2 ^ F^ ^ I transition 
of the D2 line of **^Rb, although there are also contributions from the F^ = 2 ^ Fg = 2 and 




Fig. 1 . A laser beam linearly polarized along the z axis propagates through a paraffin-coated 
vapor cell, producing ground-state alignment in the ^^Rb atoms. A constant magnetic field 
Bq and an oscillating rf magnetic field are applied within the magnetic shield that sur- 
rounds the cell. A balanced polarimeter and lock-in amplifier are used to analyze polar- 
ization rotation of the light induced by the atoms. A resonance is observed when the rf 
frequency is equal to the ground-state Larmor frequency of the atoms. 

Fg = 2^ Fe = 2> transitions. Theoretical modeling shows that signals produced on each of these 
transitions have similar line shapes, although the signal from a = 2 — >^ = 3 transition is 
of the opposite sign. In fact, the effects that we describe here are present for any transition 
with ground-state angular momentum F^ > 1. For the theoretical treatment presented here, 
we therefore consider the simplest case of a F^ = 1 — >^ F^ = transition, for which analytical 
solutions are readily obtained. The model for this system gives results in reasonable agreement 
with the experimental data, indicating that the effects are not strongly dependent on the angular 
momenta characterizing the transition. 

The atoms are subject to a z-directed field Bq — Bqz, corresponding to the Larmor frequency 
Q.L = SIJ^bBq, where jJ.B is the Bohr magneton and g is the Lande factor (we set H = 1). Lin- 
early polarized light propagating in the y direction with polarization in the z direction optically 
pumps the system and creates an aligned state. An oscillating rf magnetic field is applied in 
the X direction, Brf = ^rfxcosffltfj, corresponding to the rf Rabi frequency ilrf = gl^sBif. The 
dependence on Q,^ of the rf line shape of the optical rotation signal is studied. An analytic so- 
lution can be obtained in case in which the rf power is low enough that the rf power-broadened 
line width is much smaller than Q.l and (0^. Numerical solutions are obtained in the general 
case. 

The experiment is performed using a vapor cell with an antirelaxation coating. In such a cell, 
atoms can be optically pumped in the light beam and then exit and return to the beam after 
undergoing collisions with the cell walls, without the polarization relaxing. Thus a complete 
theoretical description must take into account the different conditions — and the different state 
of the atoms — ^inside and outside the beam. If the light power is low enough so that saturation 
effects do not occur in the beam, however, the system can be modeled by considering the 
average state of the atoms over the entire cell. This amounts to mapping the case of a coated 
cell onto that of an uncoated cell (i.e., one in which atomic polarization can be considered to 
completely relax as soon as the atoms leave the light beam). The light beam in the uncoated 
cell is taken to have an intensity equal to the average intensity over the entire cross-section of 
the coated cell, and the atomic transit rate through the beam is taken to be equal to the ground- 
state polarization relaxation rate in the coated cell. (This "effective" uncoated cell must be very 
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Fig. 2. Experimental spectra (dotted lines) and theoretical predictions (solid lines) for 
three different values of the rf field amplitude. The experimental values (obtained from 
calibrated measurements of the coil current) in terms of the Rabi frequency are (a) 
£lr{/{27t) = 3 Hz, (b) f2if/(2;r) = 63 Hz, (c) Q.,f/(2K) = 158 Hz. These values are used 
in the theory, along with the parameters natural width F = 2;r x 6.1 |is^' (chosen equal to 
the natural width of the Rb D2 transition), light power 5 |jW, ground-state relaxation rate 
7( = 2;r X 8 s^' , and bias-field Larmor frequency Q.i/{27t) = 554.75 Hz. These last three 
parameters are chosen for optimal agreement with the data — they are reasonably close to 
the experimental values. The value for light power is applied to the theoretical model using 
the method described in Sec. |2.2[ At low rf-field strengths, as in row (a), the spectra are 
Lorentzians. An additional central feature appears in the resonances when fljf exceeds y,, 
as in rows (b) and (c). Boxes in row (c) indicate regions plotted on expanded scales in Fig. 
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large, in order to account for the slow ground-state relaxation of the coated cell.) This is the 
case that we will consider. 

Another complication arising in a complete model of an atomic vapor cell is the velocity de- 
pendence of the atomic state due to Doppler broadening and coUisional velocity mixing effects. 
These effects tend to change the dependence of a signal on the light frequency. Since we hold 
the light frequency fixed in this experiment, the main effect is the introduction of an overall 
scaUng factor. Therefore, we neglect the velocity dependence in our model. 

In the Zeeman basis {\Fg = l,m = 1), \Fg = l,m = 0), \Fg = l,m = —1), |Fg = 0,m = 0)}, 
the total time-dependent Hamiltonian H of the Doppler-free system under the optical rotating- 
wave approximation is 
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where Q.r is the Rabi frequency of the optical transition induced by the linearly polarized light 
and A = 0) — cOo is the optical detuning; (O is the frequency of the light and (Oq is the frequency 
of the groimd to excited state transition in the absence of a magnetic field. For the case in 
which the ground-state relaxation rate and ilrf are both much smaller than Q.^ and (O^f, we can 
also perform the rotating-wave approximation on the rf field, in order to remove the Larmor- 
frequency time dependence from the Hamiltonian. In the rotating frame obtained using the 
unitary transformation 
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the density-matrix evolution can be written in terms of an effective Hamiltonian H' = U ^HU — 



iU ' ^ . After dropping fast-osciUating, off-resonant terms, we have 
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where Arf = ft),f — is the rf detuning. The evolution of the density matrix p (normalized so 
that Trp = 1) is described by the Liouville equation 



p = -i[H',p]-^{^,p}+A, 



(4) 



where [ ] denotes the commutator and { } the anticommutator. The relaxation of the system is 
given by the matrix 

f 7t \ 
7; 
7, 

V r+Yt ) 



(5) 



where the excited state decays spontaneously with a rate F and the ground and excited states 
relax with a rate 7, due to the exit of atoms from the Hght beam. The matrix A describes re- 
population of the ground state due to atoms entering the beam and spontaneous decay from the 
upper state, and is given by 



A = 
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where the Zeeman ground and excited sublevels are denoted as gmg and Cm^, and Pe^eo is the 
population in the excited state. 

Under the rotating-wave approximation for the rf field described above, the evolution equa- 
tions for the rotating-frame density matrix contain no explicit time dependence. We can there- 
fore set the time derivatives to zero and solve the resulting system of Unear equations for the 
steady state. Taking the case of low light power, we solve the equations to third order in the 
optical Rabi frequency, which is the lowest order at which nonlinear optical rotation signals 
appear. Using the inverse transformation U~^, we transform back to the laboratory frame to 
find the time -dependent density matrix. 

The expectation value of the medium polarization is found from the laboratory-frame density 
matrix. By substituting this value into the wave equation, we can calculate the optical-rotation 
signal measured in the laboratory frame in terms of the rotating-frame density-matrix elements. 
After multiplying by the reference signal and averaging over time, we find the in-phase and 
quadrature (out-of-phase) signals per unit length A(. of the medium to be 
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~dr ~ 4n£lR ^^(P^o^-i + P'osi ) . (8) 

where X is the transition wavelength, A' is the atomic density, and PeQg_i and PcQgi are the 
optical coherences between the excited state and the ground g-i and gi Zeeman sublevels, 
respectively. The magnitude of the optical-rotation signal is given by 



(9) 



The expressions obtained after substituting in the solution for the density matrix are compli- 
cated; to simpUfy the presentation we assimie that the Ught field is on resonance and that F is 
much greater than all other rates in the problem. This gives 

dcT ^ jVArfA^nrf(2y,^ + 8Agf - ^rf)^l 
d£ 36;rF7,0?+4A2+n2)[4(T^+A2) + n2]' 

ay""' ^ NX^^^{4f + 16A2 + ^2 )n| 
de 72;rF()^+4A2+i22)[4(^+A2)+jj2]' 



where we have neglected the contribution to optical rotation that is independent of the Ught 
power Expanding these expressions in a power series in Q.^, we obtain 



-0[£2rf]^ (14) 

These expressions describe resonances in Arf centered at A^ = 0. To lowest order in H^, they 
are proportional to the real part, imaginary part, and absolute value, respectively, of a complex 
Lorentzian. Additional features appear at higher orders, as discussed in the next section. 

When £2rf becomes of the same order as or exceeds £2/,, the rotating-wave approximation for 
the rf field is no longer valid. In this case we use the Hamiltonian H of Eq. ([T]i and proceed 
in the laboratory frame. The Liouville equation is now time dependent; periodic solutions can 
be found by expanding the density matrix in Fourier series and retaining a finite number of 
harmonics. This provides a linear system of time-independent equations that can be solved 
numerically for the Fourier coefficients. The observed optical rotation signals can then be found 
as before. 

3. Discussion 

The predictions of the density-matrix calculation described in Sec. |2.2| are compared to the ex- 
perimental data in Fig. [2] For each value of the rf field strength (characterized by the rf Rabi fre- 
quency £2rf), the in-phase and quadrature components and the magnitude of the optical-rotation 
signal as a function of the rf frequency Oif are shown. The version of the theoretical treat- 



ment valid for arbitrary rf-field strength discussed in Sec. 2.2 is used to generate the theoretical 
predictions, although for the lowest rf power, the signal is well described by the lowest-order 
terms of the expansions ([T2]i-([T4]i. Three regimes in the dependence on £2jf can be identified. 
At the lowest field strengths, £2^ < 7, the in-phase and quadrature resonances in rf frequency 
take the form of dispersive and absorptive Lorentzians of characteristic width 7, (Fig. |2^). At 
intermediate field strengths, 7 < < (we assume 7 ^ liiij), the Lorentzians broaden 
and additional narrow features are seen at the center of the resonances (Fig. |2]3), the result of 
polarization-averaging effects discussed below. 

For higher fields, £2if > |£2i|, effects due to ac Zeeman shifts and far-off-resonant fields 
are predicted to become important. We did not perform measurements in this regime, but the 
beginning of these effects can be seen in the data (Fig.|2j;). The negative-frequency component 
of therf field results in a resonance at Orf = — £2/, symmetric to the one at (Brf = ^l- In Fig |2|t') 
and (c), the off-resonant tail of this negative -frequency resonance produces an overall slope 
in the in-phase component of the positive-frequency signal. Calculations for values of Oif of 
the same order as predict higher-order resonances at odd fractions (1/3, 1/5, etc.) of £2^. 
When the experimental data presented in Fig.|2|c) are plotted on expanded scales, as in Fig. [3] 
a higher-order resonance can be observed at one third the frequency of the main resonance, in 
agreement with theoretical predictions. 

3.1. Low -field regime 

We first consider the low-rf-field regime, £irf < 7. This case was discussed in Ref |15|. An 
example of experimental data taken in this regime compared to theoretical predictions is shown 
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Fig. 3. Row (c) of Fig.|2]plotted on expanded scales. (The vertical scale is different for each 
plot.) A resonance due to high-field effects can be seen near 03^/(2k) = \^Q.i/{2Tt) ~ 185 
Hz. 



in Fig. |2|a). As described by Eqs. ([T2jl-([l4jl, the resonance observed in the in-phase compo- 
nent, quadrature component, and magnitude of the optical-rotation signal takes the form of the 
real (dispersive) part, imaginary (absorptive) part, and the magnitude of a complex Lorentzian, 
respectively. The absolute value of a complex Lorentzian has the line shape of the square root 
of the absorptive pait. Near resonance the quadrature component is the primary contributor to 
the magnitude, while farther from resonance, the in-phase component provides the main contri- 
bution. The characteristic width of the observed resonances is determined by the ground-state 
relaxation rate. 

Optical pumping by the light field removes atoms from the lUg = sublevel, leaving an 
incoherent mixture of atoms in the nig — ±1 sublevels. The atomic polarization can be illus- 
trated using the angular-momentum-probability surface (AMPS) f 26ll27ll28l . whose radius in 
a given direction is determined by the probability of measuring the maximum possible angular- 
momentum projection in that direction. This provides the quantum-mechanical analog of the 
classical angular-momentum probability distribution. The optically pumped distribution coiTe- 
sponds to atomic alignment along the z-axis with a "peanuf'-shaped probability distribution 
(Fig. |4^). (All of the AMPS shown here are obtained directly from the density-matrix calcu- 
lation. A quantity of the lowest-rank, isotropic polarization moment is subtracted from each 
figure so that the anisotropic polarization can be more clearly seen 1281.) Because the nig ~ 
sublevel has been depleted, the atomic medium transmits z-polarized light, while tending to 
absorb orthogonally polarized light — i.e., the atoms function as a polarizing filter with trans- 
mission axis along the atomic alignment axis |29 |. This linear dichroism can induce rotation of 
the Ught polarization if the transmission axis is tilted away from the light polarization axis. 

The observed signals in the low- and intermediate-field regimes can be most readily under- 
stood in terms of Larmor precession of the atomic alignment in the combined static and rf 
magnetic field. In the rotating frame, under the rotating-wave approximation for the rf field, the 
effect of the magnetic fields Bo and Brf can be described in terms of fictitious static fields Eg and 
Bjf. These fields can be determined by examining the rotating-frame Hamiltonian (|3]l. The bias 
field in the rotating frame, Bq, points along Bq, while its field strength is such that it produces a 
Larmor frequency given by the detuning of the rf field from resonance: il^ = — — —(£t^+Q.L. 
The rf field B^ becomes a static field B(.f in the xy plane with associated Larmor frequency 
= iirf/2. The direction of this field in the xy plane depends on the arbitrary phase chosen for 
the rotating frame; in our convention B[.f points along x. Thus the resultant field Bj^, = Bq + Bj^ 
lies in the xz plane. 

The atomic polarization evolves in the rotating frame under the action of Bj^f. For low field 
strengths, this evolution is Larmor precession, so that the polarization continues to correspond 
to alignment. Because all of the external fields are static in the rotating frame, the effect of 
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Fig. 4. Angular-momentum-probability surfaces illustrating the behavior of the system with 
a resonant (Ajf = 0), weak rf field, (a) The light, with electric field E, produces atomic 
alignment along the z-axis. (b) In the rotating frame, the atoms precess around the x-axis. 
Due to relaxation, a steady state is reached, (c) In the laboratory frame, the atoms precess 
around the z-axis, with the linear dichroism of the ensemble inducing oscillating optical 
rotation in the z-polarized, y-propagating light field. (Media 1 animates the precession and 
optical rotation as a function of time.) The transmission axis of the polarized ensemble, in 
this case parallel to the atomic alignment axis, is indicated with a line drawn on the surface. 

relaxation leads to a steady state for the atomic polarization. First consider the case in which the 
rf field is exactly on resonance (A^ = 0). Then B(o, is equal to B(.f and points in the x direction 
(Fig.|^), and the atomic alignment precesses in the yz plane. Under the assumption iirf ^ 7,, 
the precession frequency is much less than the relaxation rate, so that each atom precesses 
through a small angle before relaxing. The rotating-frame steady-state ensemble polarization 
thus consists of alignment at a small angle to the z-axis in the yz plane (Fig.|4|5). 

In the laboratory frame, the alignment precesses about the z-axis (Fig.|4j;). At the instant that 
the alignment is in the yz plane, it does not induce any polarization rotation in the y-propagating 
light field. On the other hand, whenever the alignment axis is tilted away from the initial light 
polarization axis (the z-axis) in the plane transverse to the light propagation direction (the xz 
plane), the atoms can induce optical rotation. Because of the precession of the alignment, the 
optical -rotation signal oscillates. The amplitude of the signal is determined by the amount of the 
alignment and the angle between the alignment axis and the z-axis (zenith angle). The phase 
of the oscillating signal is determined by the angle of the alignment axis about the z-axis in 
the rotating frame (azimuthal angle). The direction of the alignment axis in the rotating frame 
corresponds to its direction in the laboratory frame when the rf field is maximum. Only the 
component in the xz plane will induce optical rotation in y-propagating light, so an alignment 
axis in the xz plane in the rotating frame produces a signal in phase with the rf field oscillation, 
while the component in the yz plane produces a quadrature component in the signal. In the case 
of Fig.|4]the optical-rotation signal is entirely in the quadrature component. 

As the rf field is tuned away from resonance, the total field in the rotating frame begins to 
point away from the jc-axis and toward the z-axis. Precession about this field then takes the 
alignment in the rotating frame out of the yz plane (Fig. |5^). This tends to reduce the angle 
that the steady-state alignment makes with the z-axis, reducing the amplitude of the laboratory- 
frame optical rotation signal. On the other hand, because the alignment now has a component in 
the xy plane in the rotating frame, the oscillating rotation signal gains an in-phase component. 
When the rf field is tuned far enough away from resonance so that |Arf | > 7,, the precession fre- 
quency in the rotating frame becomes large enough that the atoms undergo an entire precession 
cycle before relaxing. The ensemble polarization is then averaged about the direction of the 





Fig. 5. As Fig.|4] but with an off-resonant rf field, (a) In the rotating frame, the atoms pre- 
cess around the total effective field, which lies in the xz plane. As a result, the polarization 
is no longer entirely in the yz plane, (b) Far off resonance, the effective magnetic field is 
large enough that the polarization is completely averaged about the magnetic-field direc- 
tion. (The arrow representing Bj^, is truncated, as it is too long to fit entirely on the plot.) 
Media 2 shows the polarization, along with the in-phase component of the optical-rotation 
signal, as the rf-field detuning is swept through resonant and off-resonant conditions 

total magnetic field (Fig.jsJ)). Because B[o( lies in the xz plane, the signal is now predominately 
in phase with the rf field, and the quadrature component is strongly suppressed. As the detuning 
becomes large, the average polarization points more and more along the z-axis, and the signal 
amplitude drops to zero. 

The preceding description is seen to correspond to the signals shown in Fig.|2|a) — the dis- 
persive (in-phase) and absorptive (quadrature) components and the magnitude of a Lorentzian 
with characteristic width 7, — as well as to the lowest order terms of Eqs. ([T2]i-([T4]i. 

3.2. Intermediate-field regime 

When the rf field is large enough that Jt < Clrf < the rotating-frame precession frequency is 
high even at zero detuning. This causes averaging of the atomic polarization about the magnetic- 
field axis. For — this is the jc-axis: polarization transverse to the jc-axis is averaged out. 
However, the polarization along the x-axis is preserved, so that the x-axis becomes the preferred 
axis for the polarization (Fig. |6^). The "doughnut"-shaped probability distribution seen in Fig. 
[6|a) is obtained from the initially pumped "peanut"-shaped distribution (Fig.|4^) when copies 
of the peanut distribution rotated by arbitrary angles about the x-axis are averaged together. 
Another way to explain the doughnut shape is to transform to the basis in which the quantization 
axis is along x. In this basis, the excitation light is a polaiized, so that it pumps atoms out of the 
bright state consisting of a superposition of the m^ = ± 1 sublevels, and leaves them in the dark 
state made up of the opposite superposition, as well as in the mg = sublevel. However, due 
to the precession induced by the x directed magnetic field, atoms oscillate between the bright 
and dark superpositions, so that the pump light removes atoms from the m^ = ±1 sublevels 
incoherently. The atoms are then left in the m^, = sublevel, i.e., the atoms have no angular- 
momentum projection on the x-axis and are symmetric about the x-axis, as seen in Fig.|6|a). 

Polarization along the x-axis remains in the xy plane as it precesses around the z-axis in the 
laboratory frame and so does not induce any optical rotation. However, when the rf field is tuned 
slightly away from resonance, Bj^, points away from the x-axis, and so also does the averaged 
atomic polarization (Fig.jSJ)). The polarization then causes optical rotation. As the polarization 
in the rotating frame is in the xz plane, the signal in the laboratory frame is in phase with the rf 
field. (The plots in Fig. |6] are shown for the instant at which the rotating frame coincides with 



z z 




Fig. 6. Rotating-frame AMPS for intermediate rf field strength, in the ideal case in which 
7r fljf <C Optical rotation is indicated for the instant in which the rotating frame 
coincides with the laboratory frame — this means that the generation of the in-phase com- 
ponent of the signal is shown, (a) When the rf field is on resonance, the atomic polariza- 
tion (created along the z-axis) is averaged about the x-axis. The resulting polarization is a 
"doughnut" aligned along the x-axis. The transmission axis, perpendicular to the alignment 
axis, is marked with a line on the surface, (b) For small detunings, the alignment axis fol- 
lows the effective magnetic field direction, (c) When the effective magnetic field is at the 
magic angle 9,„ = arccos(l/\/3) to the light polarization direction, the atomic polariza- 
tion is completely averaged out due to precession, (d) For larger detunings, the polarization 
regains its original "peanut" shape, and the transmission axis is along the alignment axis. 
Media 3 shows the atomic polarization and in-phase optical rotation as the rf-field detuning 
is swept through resonance. 



Fig. 7. Cross section of the AMPS for pure alignment along the z-axis, described by 
Y2fl{0,<j)) °= 3cos^ 9—1. Positive function values are shown as solid lines, negative func- 
tion values as dashed lines. The magic angle 6,„ = arccos( 1 / \/3) is indicated. 

the laboratory frame, so that the in-phase component of the optical rotation is shown.) 

To understand the generation of optical rotation in this case, it is important to note the 
effect of the different character of the ensemble polarization. In the low-field case, the 
angular-momentum probability distribution has maxima along the ensemble polarization axis 
(a peanut), while in the current case, the distribution has minima (a doughnut). (This result de- 
pends on the particular type of transition considered — for a Fg Fe=Fg transition the situation 
is reversed.) We can think of the doughnut shape as an unpolarized distribution (sphere) with 
a peanut shape removed. From this viewpoint, it is reasonable that this "negative polarization" 
produces rotation of the opposite sign. More concretely, we can note that, as described above, 
in a doughnut distribution atoms are concentrated in the nig = sublevel with the quantiza- 
tion axis along the alignment axis. The state then preferentially absorbs light that is polarized 
along the alignment axis. In the analogy with a polarizing filter, the transmission axis of the 
doughnut-shaped probability distribution is transverse to, rather than along, the alignment axis. 
This explains the sign of the rotation shown in Fig.[6|b) when B((„ points away from the x-axis. 
(The transmission axis of each polarization state is marked with a line on the surfaces plotted 
inFig.|6]) 

As the rf field is tuned farther from resonance, Bj^j points farther away from the x-axis, 
bringing the averaged alignment axis with it. This larger angle produces a larger optical rota- 
tion signal. However, another trend eventually takes over: because of the shape of the initially 
pumped polarization distribution, the amount of polarization that lies along the magnetic-field 
direction decreases. As a result, the amount of averaged polarization is reduced, tending to 
reduce the signal. In order to analyze this, we can plot just the aligned part of the initially 
pumped density matrix, neglecting the isotropic part that is included in Fig. |4|a). The surface 
corresponding to just the rank K — 2, q = Q polarization moment is described by the spherical 
harmonic Y2X){Q,i^) °= 3cos^ 0—1, plotted in cross section in Fig. [t] Negative values of the 
function are indicated by dashed lines. There are maximum positive values along z (0 = 0), and 
maximum negative values in the xy plane (0 = n/2). As 9 moves away from either of these 
values, the magnitude of the polarization is reduced. At a particular angle 0,„ = arccos(l /\/3), 
analogous to the magic angle observed in nuclear magnetic resonance experiments, the polar- 
ization moment goes to zero. This means that if Bj^, is at this angle to the z-axis, the averaged 
polarization completely cancels, sending the optical rotation signal to zero (Fig.|6};). This con- 
dition corresponds to the additional zero crossings seen in the in-phase component of Fig.|2|b,c) 
above and below the center of the resonance. The direction of the effective magnetic field in the 
rotating frame depends on the rf field strength and the detuning from resonance — as the field 
strength is increased, the detuning required to achieve the magic angle also increases. 



When the detuning is even larger, there is once again residual polarization after averaging 
about the magnetic-field direction. However, now the magnetic-field direction is close enough 
to the initial alignment axis that the polarization resulting from averaging resembles the peanut 
shape of the initially pumped polarization (Fig.|6}l). Thus, the optical-rotation signal in this case 
is opposite in sign to that for small detuning and has the same sign as that for the low-power 
case. 

As the detuning continues to increase, B't^t and the averaged atomic alignment point more 
toward the z-axis, reducing the optical-rotation signal. 

The preceding discussion describes the in-phase signal shown in Fig.|2|b) and (c): a power- 
broadened Lorentzian with a narrower feature of the opposite sign in the center. If the discussion 
is strictly interpreted, there should be no quadrature signal in this regime, as the polarization 
in the rotating frame is always in the xz plane. Figure |2] does display (strongly suppressed) 
quadrature signals, which are a remnant of the low-field regime. 

3.3. High-field regime 

As iljf becomes of the same order as £2^ or exceeds it, various higher-order effects appear in the 
data and the full theory that cannot be described under the rotating-wave approximation for the 
rf field. In particular, a resonance near £2^/3 is seen, as shown in Fig. [3] This can be explained 
as due to ac Zeeman shifts, which produce evenly spaced sidebands that result in resonances at 
odd subharmonics of the lowest-order resonance. The additional features can also be interpreted 
as arising from higher-order resonances between the Larmor precession and the rf frequency, 
similar to those seen in nonlinear magneto-optical rotation with frequency-modulated light f30l . 
As the rf field strength increases, many additional resonances are predicted by the theory. 

4. Conclusion 

We have conducted a detailed experimental and theoretical investigation of magneto-optical 
rotation in the presence of a strong radio-frequency field. A model has been developed that 
allows both the reproduction of the nontrivial line shapes that appear as the rf field strength is 
increased, and also the qualitative understanding of the atomic polarization dynamics respon- 
sible for the signals. The results of this study will be useful for the analysis of the behavior 
of atomic magnetometers exposed to strong radio-frequency fields. The results may also have 
applications in the design of efficient methods for the preparation of specific polarization states 
of atoms, for example, atomic states with large alignment but no orientation, as are used in the 
workof Ref. lfT9l . 
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